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Abstract 

In the Laguerre ensemble of N x N hermitian matrices, it is of interest both 
theoretically and for applications to quantum transport problems to compute 
the variance of a linear statistic, denoted var^r /, as N — > oo. Furthermore, 
this statistic often contains an additional parameter a for which the limit a — * 
oo is most interesting and most difficult to compute numerically. We derive 
exact expressions for both limjv-»oo var^r / and lima^oo liniAr^oo var^r /■ 

I. INTRODUCTION AND STATEMENT OF MAIN RESULT 

In the random matrix theory of quantum transport (see 0,0,0 and references therein) 
the following quantity is fundamental: 1 



var N f : = jf f 2 K N (ji, fi) dfi 

-ff'H'ffl*' ** (L1) 

where K N (fi,fi') is the Laguerre kernel; that is 

JV-l 

K N {fi,fi') = ^(/^V) 

3=0 
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lr The notational choice f(y/x) rather than f(x) in (1.1) will be convenient later. It also agrees 
with the convention of Stone, et. al. [Hi. 
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and {(pj(x)} is the sequence of functions obtained by orthonormalizing the sequence 

I J j=0 

over (0, oo) (here v > —1). In particular, one is interested in 

var / := lim va.r N f (1.2) 

in the limit a — > oo. In applications various choices are made for /, but we need assume 
here only that / is smooth and sufficiently decreasing at infinity to make the integrals 
well-defined. 

In the random matrix model of disordered conductors, the quantity var / is related 
(via the two-probe Landauer formula) to the fluctuations of the conductance and the limit 
a — > oo is the high-density (or metallic) regime. A lucid account can be found in the 
review article by Stone, Mello, Muttalib, and Pichard || to which we refer the reader for 
further details and references. However, these authors did not evaluate varjv/ in the limits 
of interest; namely iV — > oo followed by a — > oo. It is the purpose of this paper to evaluate 
these limits. We will see that the result agrees with the prediction of Beenakker |E]|| who 
gave a heuristic argument for this limit. 

By a change of variables we write ( |1 . 1|) in the more suggestive form: 



- 1: i: ^W){^Mw^)P^ <-> 

From asymptotic formulas for generalized Laguerre polynomials (see, e.g., 10.15.2 in ||) it 
follows that (§|7| 

K(x,y) : = lim Km{——,-^—) 

Ju(\fx) y/yJ' u (y/y) - y/xJ' v (y/x) Jy{^/y) , 

2(x-y) 1 ' ] 

where J v {z) is the Bessel function of order v. (The limit is uniform in x and y for < x, y < 
L < oo and all L.) We call K(x, y) the "Bessel kernel." (This kernel also arises in scaling 
the Jacobi ensemble of random matrices at either edge ±1.) Using this in ( |1.3| ) we obtain 

var f = [ f 2 ( — — ) K(x, x) dx 
Jo \ a J 

- f f f (^f\ f (^pj K 2 (x, y) dxdy (1.5) 

where K is the Bessel kernel. 

The problem is reduced to evaluating ( |1.5|) in the limit a — > oo. We will show that 



1 r 00 
lim var / = — - / 

a ^°° ' n 2 J-oc 



f(2iy) 



y tanh(-7n/) dy (1.6) 
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where / is the Mellin transform of /, i.e. 

f{z) = / x z - l f(x)dx. 
Jo 

This agrees with the result of Beenakker once one notes that his f(x) is our f(y/x). 
For numerous applications of ( |1.6| ) we refer the reader to Beenakker M. 



II. THE LIMIT a -» oo 
A. Use of Hankel Transform 

It is convenient to define the kernel 

„2 „,2\ 



L(x,y) : = 2K{x\y< 

tJ v (xt)J v (yt)dt (2.1) 



where K(x, y) is the Bessel kernel. (A simple proof of the second equality can be found in 
0.) Then var / can be written 

var / = J™ xf 2 (^JL{x,x)dx-h (2.2) 

where 

h 



I I xyf ( - ) / (-) L 2 (x } y) dxdy 
Jo Jo \aj \aj 



i ,oo / ri / roo ^ x ^ ^ ^ t dAyf(^)tJ v (ty) dydt 

o Jo \Jo \Jo a ) ) a 

i ,00/ r/ r j v ( te )j^^) dA'j^) ^/(^tj^y) ^ 

o Jo \Jo \Jo ) a J a 

(2.3) 

where we used ( [2.1|) to deduce the middle equality and we made the change of variables 
x/a — > x and at' — > t' to obtain the last equality. 
We now recall the Hankel inversion formula: 

u\ / xg{x)J u {xu)dx)J u {u£ s )du = g{£ t ) 



which hold for y/xg(x) continuous and absolutely integrable on the positive real line and 
v > — |. First writing the t'-integration in (|2.3|) as the integral from (0, oo) minus the 
integral from (a, oo) and then employing the Hankel inversion formula (with the choice 
g(x) = f{x)J u {atx) ) on the part containing the t'-integration from (0, oo), we see that this 
part exactly cancels the single integral appearing in the expression (|2.2|) for var /. Thus we 
are left with 
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rl roo roo roo y jjy 

var f = dt dy dt' dxxytt' f(x)f(—)J v (atx)J u ( — )J u (t'x)J u (ty) 
Jo Jo Ja Jo a a 

rl roo roo roo 

= a A I dt I dy I ds dxxytsf{x)f{y)J v {atx)J u {asx)J v {asy)J v {aty). 
Jo Jo Ji Jo 

(2.4) 

We remark that the Hankel transform plays the analogous role for the Bessel kernel that 
the Fourier transform plays for the sine kernel 

1 sin 7r(x — y) 
7T x — y 

in the Gaussian Unitary Ensemble. 



B. Residue Calculation 



Introducing the (inverse) Mellin transform 

I rc+ioo 



rc+too „ 

f(x) = — f(z)x- z dz (c>0) 

Z1X% J c—ioo 



into ( |2.4| ) and interchanging the orders of integration we see that 

I rc+ioo A ] rc+ioo 

var / = a 4 — / dz x f{z\)— / dz 2 f(z 2 ) 

Z7TI J c— ioo Z7TI J c—ioo 

rl roo roo roo 

x / dtt / ds s / dxx~ Zl+1 J l/ (atx)J 1/ (asx) / dyy~ z ' 2+1 J u (aty)J u (asy). 
Jo Ji Jo Jo 

The x and y integrations can be performed using (6.5762) in ||; namely, 



roo 

I x~ x J u {ax)J v {bx)dx 
Jo 



2 



2 A (a + &)^- A+1 r(i + z/)r 



l+A 
2 



/ 1 - A 1 4ab \ 

xF ^ + _,„ + _ ;2 „ + 1; __j 

where F(a, b; c; z) is the hypergeometric function, a, b > 0, 29ft(z/) + 1 > 3?(A) > — 1. 
In the resulting integral we make the following change of variables: 

Ats 

u = u i ^2 ' v = t + s 
{t + s) 2 

which has Jacobian 

11 

J(u, v) 



In the v u-plane we are now integrating over the region in the first quadrant bounded above 
by the curve 



4 



4(v- 1) 



v > 1 



and bounded below by the ray [1, oo] on the v-axis. The v- integration may now be trivially 
done with the result that 



var / 



42«,-ip2( I/ + x) 2m 
1 



V (v + 1 - z x /2) 



r (^/2) 



x ■ 



2tT2 



c+io ° , // \ /(H1-W2) 



r (V2) 



X- 



Zl + Z 2 Vo 



-1/2 



u 



Zl+Z 2 



(l-vT 



it 



21+22 



where 



i^fou) :=F(i/ + 1--,i/ 



^;2z/ + l;w 



We now use zi and z = z\ + -22 as integration variables so that 

1 1 



var / 



1 



' +- °°^/> 1 ) r( ?, 1 7! l/2) 

c— ioo 1 \Z\ 1 1) 



X 



27U J2c-ioo 
•1 



dz/(z - zi)a 



T((z- Zl )/2) 



x 



■it 







-1/2 



1 + v 7 !^ 



it 



xF u (zi,u) F u (z - zi,u) . 

Observe that the a dependence of var / resides solely in the term a z in the above integral. To 
compute a — > 00 this suggests we should first deform the contour into the left-half z-plane. 
The lim^oc var / will be determined by the residue of the pole at z = 0. 

To calculate this residue (which is a function of z\) we must know the principal part 
of the Laurent expansion (in z) of the integral involving the ^-integration. The divergence 
of this integral as z — > is determined by the behavior of the integrand in the vicinity of 
u = 1. This behavior near u = 1 is straightforward to compute since it is known [Q that 



F(a, b; c, u) 



T{c)T{a + b-c) 

mm 



\ c—a—b 

u) as u 



1. 



Thus lim^oo var / is expressed as a single integral over the variable z\ . If we now make 
use of the T-function identities [HI 



T(z)T{-z 

r(i/2 + *)r(i/2-z 



7T 



zsmnz 

7T 



COS TTZ 

F(2v + 1) = 2 2u tt- 1/2 T{u + 1/2)I> + 1) 
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we obtain 

Jim var / = -_ f(z 1 )f(-z 1 )-t B! a(-zOdz 1 . 
We now deform the contour to the imaginary axis (and send y — > 2y) to obtain (|l 
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